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We report on a measurement of the characteristic temperature of an ultracold rubidium ion source, in which 
a cloud of laser-cooled atoms is converted to ions by photo-ionization. Extracted ion pulses are focused on a 
detector with a pulsed-field technique. The resulting experimental spot sizes are compared to particle-tracking 
simulations, from which an effective source temperature T = (3 ± 2) mK and the corresponding transversal 
reduced emittance e r = 1.4 x 10~ 8 m rad y/eV are determined. Space charge effects that may affect the 
measurement are also discussed. 
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PACS numbers: 29.27.Bd, 32.80.Fb, 41.75.Ak, 67.85.-d 



I. INTRODUCTION 

Focused Ion Beams (FIBs) are widely utilized in the 
semiconductor industry and in nanosciencei^. They are 
used successfully for high precision milling or deposition 
and on the other hand, similarly to scanning electron 
microscopy (SEM), for high resolution microscopy^. To 
keep up with reduction in sizes in the semiconductor in- 
dustry it is necessary to reduce the smallest focusable 
spot size^£. The properties of the source have a key role 
in what can be achieved. A typical FIB using the cur- 
rent industry-standard Liquid Metal Ion source (LMIS) 
can reach a high brightness (10 6 A m~ 2 sr _1 eV _1 ) and 
can deliver 10 pA current in 10 nm spot size^ (note that 
higher currents can be reached if the beam is less fo- 
cused). Brightness is the current per unit area and solid 
angle, normalized by the beam energy and it is a key 
property for a source. A new kind of ion source, the 
ultracold ion source (UCIS) has been proposed as an al- 
ternative for the LMIS£~— . The UCIS is based on creating 
very cold ion beams by near-threshold photo-ionization 
of a laser-cooled and trapped atomic gas, with a source 
temperature T less than 1 mK. The UCIS has the poten- 
tial of producing ion beams with a brightness of 10 5 A 
m -2 sr _1 eV _1 and a current up to 100 pA, according 
to simulations^. Its major advantage is an energy spread 
two orders of magnitude lower than the LMIS (down to 
10 meV) as demonstrated by Reijnders et a/3. Lower 
energy spread may lead to a smaller achievable spot size 
by reducing the contribution of chromatic aberration!. 
This work aims to measure the source temperature of a 
rubidium UCIS, which is an important physical quantity 
related to the ability to focus the beam. Here we show 
how an effective source temperature can be extracted 
from measurements of the minimum spot size achieved 
in focusing the beam and we demonstrate a new way to 



focus ion beams. The ultra low temperature of the source 
permits collimated beams to be created at low energy 
(down to a few eV— ), which allows using time-dependent 
fields for accelerating and focusing. The duration and 
the shape of the accelerating electric field pulse can be 
tuned so that a variable focusing lens is created. An ef- 
fective source temperature is then extracted from waist 
scans varying the focal strength of the time-dependent 
lens. 



II. PRINCIPLE OF THE MEASUREMENT 

In the experiment, ion bunches, extracted from the 
source with longitudinal energy U, are focused on a de- 
tector with a variable focal length /. With a simple first 
order optical model it is possible to find an upper limit for 
the source temperature from the final root-mean-square 
(rms) spot radius a Xf (see Fig. [1] for a schematic draw- 
ing), assuming that the emittance is conserved. 



+ 



Ox [ 




I f 


L"l 











a ' Electronic mail: E.J.D.Vredenbregt@tue.nl 



FIG. 1. Principle of the source temperature measurement. 
Ion bunches are focused on a detector with a time-dependent 
lens with focal length /. The quantity a x i is the initial rms 
angular spread, a Xi is the initial rms radius of the source and 
a Xf the rms radius of the final spot on the detector. 
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The size of the ion bunch at the detector position is 
completely determined by the initial rms angular spread 
<r x ' when the image is in focus, 



&X f — j ®x' x • 



(1) 



The rms value of the angular spread is related to the 
source temperature T by 



2U ' 



(2) 



where kb is Boltzmann's constant and U is the energy of 
the bunch. Reversing this equation, the source tempera- 
ture can be defined given a measurement of a Xf as 
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(3) 



Any emittance growth, e.g., due to distortions in the 
imaging system, will affect the extracted value of T, 
which thus becomes an effective source temperature. 
Knowing T and thus o x \ , we can derive the rms nor- 
malized emittance e r of the source when, in addition, the 
source radius a Xi is known. The normalized emittance 
is a measure of the phase space occupied by the beam 
multiplied by the square root of the beam energy and is 
given by"i£ 



Indeed in Eq. (UJ) the dependence on the energy U dis- 
appears indicating that the normalized emittance is a 
property of the source. 



III. EXPERIMENTAL SETUP 

The UCIS is based on the technique of laser cooling and 
trapping of neutral atom o 11 ' 12 . We start with a magneto- 
optical trap (MOT) of rubidium atoms. Three orthog- 
onal pairs of counter propagating 780 nm laser beams 
(trapping laser beams) are used to Doppler cool a 85 Rb 
atomic cloud and a quadrupole magnetic field is added 
to trap the atoms. Typically 10 s Rb-atoms are trapped 
in a volume of 16 mm 3 rms at an expected source tem- 
perature of T = 143 fiK n , which corresponds to 9 neV 
of average kinetic energy of the atoms. The MOT is sur- 
rounded by an accelerator—. 

The accelerator has a cylindrically symmetric structure 
and it is placed in a vacuum chamber where the rubid- 
ium pressure is 10~ 9 mbar. See Fig. [5] for a schematic 
drawing of the experimental setup. The trapping laser 
beams enter through openings present in the structure. 
The atomic cloud is trapped in the middle of the accel- 
erator at the intersection of the trapping laser beams, 10 
mm away from the accelerator's exit, which consists of 
a circular hole with a diameter of 20 mm. The electric 



field strength at the starting point zn is 0.37 kV/cm per 
kV input voltage V a and U = eV a /2.05, where e is the el- 
ementary charge. For details of the accelerator structure 
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FIG. 2. Schematic view of the experimental setup (not to 
scale) in the x-z plane. The ionization laser beam and the 
excitation laser beam select a portion of the Rb atomic cloud 
trapped inside the accelerator. The bunched ions are acceler- 
ated using a time-dependent potential V a (t) and after a flight 
distance of 1.51 m they reach a MCP detector with a phosphor 
screen. The images are captured by a CCD camera. 

The ionization mechanism is a 2-step process. In order 
to avoid spherical aberration in the focusing process, it 
is necessary to work with a smaller volume than the to- 
tal MOT. To achieve this, the trapping laser beams are 
turned off for 25 fis. During this time, a laser beam (ex- 
citation laser beam) with the same wavelength (780 nm) 
is focused horizontally along the z-direction to a rms ra- 
dius (Jexcit = 54 /im. The excitation laser beam excites 
the Rb cloud to the 5p-level for 2 /is. Coincidentally, a 
479 nm ionization laser beam is sent in vertically along 
the x-direction for 400 ns. The ionization laser beam also 
has an rms radius <Ji on i Z = 54 /im. Fig. [3] shows an ex- 
perimental example of the timing sequence used for the 
ionization of a portion of the Rb atomic cloud trapped 
in the MOT. In this way, only the atoms at the intersec- 
tion of the 2 laser beams will be ionized. The shape and 
position of the ionized cloud can be changed by changing 
the size of the lasers and their position. The position and 
dimension of the lasers are controlled by 2 CCD cameras, 
virtually positioned in the center of the MOT, which is 
at zq on the beam axis. The initial ionization volume 
is not spherical: in the x-direction the rms radius a x is 
determined by the intersection of two laser beams and 
hence 
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so that a x is a factor y/2 smaller than the size of the two 
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FIG. 3. Typical example of the timings used in the experiment 
to ionize a portion of the whole atomic cloud. The trapping 
laser beams are turned off for about 25 /is and during this 
time the excitation and ionization laser beams are turned on 
in coincidence. The excitation laser beam is turned on for 2 
fj,s and the ionization laser beam for 400 ns. 



laser beams. 

The initial ionization volume size has been experimen- 
tally and numerically determined as follows. Ion bunches 
arc created in a DC electric field and arc consequentially 
accelerated towards the detector (which is described be- 
low). With a DC accelerating field, the exit of the acceler- 
ator forms an aperture lens with a negative focal length 
fo of 33 mm ("exit kick" effect), which is independent 
of the acceleration voltage^. This leads to a magnified 
image of the source on the detector with a known mag- 
nification of 46. From the final spot size, it is therefore 
possible to determine the source size. The initial angu- 
lar spread has a negligible influence (less than 1%) on 
the final spot size for U > 200 eV and this measurement 
was performed at U = 3 keV. The initial volume is as- 
sumed to be Gaussian distributed in all 3 directions from 
the fact that the profile of the lasers is Gaussian as well. 
Analysis of the recorded detector images give a x — (38 
± 2) /im and a y = (54 ± 3) /im. In the z-direction, a z — 
54 /im because it equals the rms radius of the ionization 
laser beam <Ji on iz ■ 

A double multichannel plate (MPC) detector with phos- 
phor screen is located at 1.51 m from the ion's starting 
position zq. The diameter of the MCP is 40 mm. A 16 
bit, cooled CCD camera is used to image the phosphor 
screen through a lens placed in between with a magnifica- 
tion of 0.33. The spatial distribution of the ion bunches 
can be extracted from the CCD images. The resolution 
of the detector has been determined experimentally by 
placing two pinholes, with diameters of 25 /im and 50 
/zm, downstream in front of the MCP. Ion bunches with 



U = 3 keV passed through the pinhole. The rms radius 
a a of an aperture with diameter D is equal to D/4. The 
rms spot radius <7<j e t measured at the detector was sub- 
stantially larger than (7 a . clS cL result of the resolution of 
the detector. Multiple experiments, with both the pin- 
holes, have been performed in order to obtain sufficient 
statistical accuracy. Using 



Vdet 



(6) 



the rms resolution of the detector 5 is found to be (95 ±4) 
/im. In what follows, S is quadratically subtracted from 
the measured spot in order to obtain a XJ = y /(T 'iet ~ ■ 



IV. FOCUSING WITH TIME-DEPENDENT ELECTRIC 
FIELDS 

Measuring the effective temperature of the source re- 
quires focusing the ion bunches onto the detector. In 
fact, as shown in Section UU at the beam waist the final 
spot size is completely determined by the initial angular 
transverse spread, which depends on the source temper- 
ature. To focus the ions, a positive lens is required, while 
the aperture lens formed by the accelerator is negative. 
However, the sign can be reversed using time-dependent 
accelerating fields. Due to the low temperature of the 
ions, it is possible to form well-collimated ion beams at 
very low energies, U ~ 10 — 50 eV. Such ions travel for 
few microseconds in the accelerator before exiting it. For 
instance, in the case of U — 30 eV, the traveling time is 
in the order of 3 /is. This time is long enough to switch 
the accelerating field while the ions are still in the accel- 
erator. 

When a constant positive voltage is applied to the an- 
ode the ions will experience a negative lens effect. The 
on-axis longitudinal electric field E z (z) follows approxi- 
mately a Gaussian centered on z , see Fig. H](b). Fig. @] 
shows the potential applied to the anode V a (t), the lon- 
gitudinal electric field component along the beam axis 
E z (z) and on-axis radial electric field component E r (z) 
for r 7^ 0, in case of a static (left hand side) or bipolar 
field (right hand side). The cylindrical symmetry of the 
accelerator makes it possible to write the electric field in 
the accelerator as an expansion of the electric field on the 
symmetry axis. In a first order approximation, the radial 
electric field component E r (r, z) is given by 



E r (r,z) 



r dE z 
2 dz 



(7) 



where r is the radial position. In the case of a static elec- 
tric field, the focal strength of the accelerator is described 
by 
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where fo is the time-independent focal length, and the 
static field e z (z) = E z (z)/<f>o, with the potential given by 
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FIG. 5. Example of a typical bipolar anode voltage pulse V a 

- 125 V is the positive 



used in the experiment. Here, V p , 
voltage with duration t v — 2 /is and V„. 



-90 V is the 



negative voltage with duration t„ = 4 /is. The pulse is created 
with a programmable waveform generator, amplified 50 times 
and applied to the accelerator. 



FIG. 4. The potential applied to the anode V a , the longitu- 
dinal electric field component E z and the radial electric field 
component E r (for r ^ 0), in case of a static (left hand side) 
or bipolar field (right hand side) . The quantities r p and z s are 
the length respectively in time and space of the positive part 
of the pulse. In panels (d), (e) and (f), the thin line indicates 
the static case. The y-axis is in arbitrary unit in all the six 
plots. 



0o = E z (z')dz' : as suggested inii. Because the inte- 
gral of the radial electric field experienced by an ion over 
time is positive (see Fig. |4] (c) in dark), a defocussing 
effect results. 

By turning off the electric field before the ions leave the 
accelerator one can cancel the exit kick effect. Moreover, 
with the use of a bipolar pulse, the exit kick effect can 
be reversed and a focusing lens created, see Fig. 0] (f). 
In this case, the net integral of the radial electric field 
component will then result to be negative. Typically the 
electric fields can be changed in less than 100 ns. Bipolar 
pulses are created with a programmable waveform gen- 
erator and amplified 50 times. An example of a bipolar 
pulse V a (i) used in the experiment can be seen in Fig. [5] 
Here we can define the positive voltage V pos with dura- 
tion t p and the negative voltage V neg with a duration r n 
long enough in order that the ions have already left the 
accelerator when V a goes to zero. The accelerating pulse 
is turned on within 100 ns after the ionization laser is off. 
The ion bunch moves downstream in the direction of the 
detector (see Fig. [2]). 

Now a "waist scan" can be performed, i.e. the final spot 
size is measured depending on one of these parameters. 
What is actually varied then is the focal strength 1/ ft, 



which in the case of a bipolar pulse, can be expressed as 
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Here, z s is the position to which the center of the ion 
bunch has moved when the field is switched from positive 
to negative (at time r p ). Obviously, the focal strength 
depends on V pos , V neg and r p , not only directly but also 
through their influence on z s . In the experiment we only 



vary V ne g, for a fixed Vp OS 
given by 



U 



and t p . The bunch energy is 



E z dz. 



(10) 
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It is of interest to note that U varies during a waist scan 
due to the different values of V neg ■ 

For further details about focusing with time dependent 
fields see Reijnders et al~. 



V. ANALYSIS 

Images collected from the CCD camera are fitted to a 
2-dimensional Gaussian. Due to the fact that the initial 
volume is not spherical (see Section IITl]) . this will result 
in a non-circular spot on the images. The longer size 
is by definition in the y-direction and the short size is 
in the x-direction. Only the x-direction is considered in 
this paper because along y-direction the bunch suffered 
from aberrations, which we ascribe to the long beam line 
and low energy of the ions, which makes them particu- 
larly sensitive to small external fields. In the x-direction, 
residual effects cannot be completely excluded, in which 
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case the extracted temperature will effectively include 
such distortions, which will generally increase the source 
emittance. 

Particle tracking simulations are performed with the 
General Particle Tracer code (GPT)^£ in order to repro- 
duce the measurements and fit the data. The electric field 
inside the accelerator has been calculated with the Su- 
perfish Poisson solver—. The initial particle distribution 
is a 3D Gaussian with dimensions a x , a y and 172; as listed 
in Section \UJ\ centered at position zq. The initial veloc- 
ities follow a Boltzmann distribution. A time-dependent 
V a (t) is applied depending on parameters varied exper- 
imentally, i.e. V neg . The rms radius of the simulated 
images are calculated and compared to the experimen- 
tal data. Minimization of the reduced y 2 is performed 
when comparing the measurement data with GPT sim- 
ulations in order to extract a parameter such as the ef- 
fective source temperature. According toi£, the accuracy 
dp of the parameter that is optimized by minimizing the 
X 2 can be calculated from the dependence of x 2 on that 
parameter in the region of the minimum as 
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FIG. 6. Spot radius a x , versus negative pulse amplitude Vneg 
for 4 different V pos (indicated in the figure) . The bunch energy 
U at the minimum of the waist scan (which depends on V poa 
and V neg ) is also indicated in the label and varies from 14 eV 
to 34 eV. 
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VI. RESULTS AND DISCUSSION 

The measurements presented here consist of waist 
scans where the negative voltage of the time-dependent 
pulse V neg has been varied. Because of the change of 
V neg , the focal strength of the time-dependent lens varies, 
as shown in Eq. ([9]). The other parameters are fixed: 
t p = 2 (is, r n = 4 [is (see Fig. [5} and V pos = 125 V (ex- 
cept in the measurement in Fig. [51 where this parameter 
is also varied). 

The repetition rate during all the experiments is 20 kHz 
and the exposure time of the CCD is 2 s. Thus every 
image is the sum of 40000 bunches. 
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FIG. 7. Final spot radius squared cjL versus the reciprocal 
of the bunch energy 1/U . The thick line is a linear fit, whose 
slope is proportional to the effective source temperature T. 



A. Temperature determination 

Fig. |5] presents four waist scans at four different V pos . 
For each V pos , the minimum spot radius occurs at differ- 
ent Vneg, and according to Eq. (fl"0|) also at a different 
beam energy U. The plot shows that the minimum spot 
radius <j xj increases when U is lowered: lower energy 
means longer time of flight and this results in a larger 
final spot radius. This also indicates that the measure- 
ment is indeed sensitive to the effective source tempera- 
ture. To find T, first the minimum spot radius is calcu- 
lated by fitting the bottom part of a waist with a second 
order polynomial. The extracted minimum spot radius 
squared a 2 is plotted versus U~ x in Fig. [71 According 
to Eq. ©, T can be extracted from a linear fit. The ef- 
fective source temperature is found to be T = (4.9 ± 0.3) 



mK. This temperature, while it is extremely low, is still 
rather high compared to the expected source tempera- 
ture (To = 143 j«K), even when To is corrected for the 
fact that the ionization laser beam was turned 0.6 nm 
above the ionization threshold. In fact, this would make 
T,5 = 390 (jX. 

To minimize any effect of space-charge forces on the ef- 
fective source temperature, we analyzed in detail a waist 
scan obtained at en even smaller charge of about 0.022 
ions per bunch (on average). The data is shown in Fig. 
ID The effective source temperature can be immediately 
estimated from Eq. ([3]), as done in the previous para- 
graph. Using U — 32 eV (the bunch energy at the min- 
imum position) and a Xf = (73 ± 4) /xm, we find that 
the resulting effective source temperature T = (1.8 ±0.2) 
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FIG. 8. Waist scan performed measuring the final spot size 
a Xf versus the negative voltage of the time-dependent electric 
field Vneg- The experimental data is scattered with the error 
bars. The solid and the dashed lines are two GPT fits which 
do and do not include space charge forces respectively. 

mK. This first order approximation is confirmed by a fit- 
ting procedure, where the effective source temperature 
can be extracted more precisely by fitting the behavior 
of a XJ versus V neg (Fig. [8]) with particle tracking sim- 
ulations. Waist scans are simulated for several source 
temperatures and the ^-minimization procedure from 
Section [V] is applied. Fig. [8] shows a GPT simulation 
with parameters that best minimize the \ 2 f° r different 
T (solid line). The extracted value is T = (3 ± 2) mK. 
The value of x 2 — 1-5 indicates that the fit is good but 
could be further improved. 

B. Influence of space charge 

The measurements of Fig. [5] were taken using far less 
than one ion per bunch and high repetition rate. There- 
fore, space charge forces should not be important in this 
case; but it is nevertheless interesting to investigate at 
which level Coulomb interactions can still play a role. 
The number of ions which are ionized by a laser beam 
follows Poissonian statistics, so it is possible, even when 
the expected average number of ions per bunch fj, is less 
than one, to have some bunches with 2 or more ions. The 
Poisson distribution is given by p n (/j,ri) = (/i"e"'' /n\), 
where n indicates the number of ions per bunch. As an 
example, when the expected number of ions per bunch is 
0.022, the probabilities to obtain one, two or three ions 
per bunch are respectively 1.96%, 0.02% and 0.0001%. 
The probability for zero ions is the highest (about 98%) 
since only one ionization laser pulse every 50 ionizes one 
or more atoms. GPT simulations for n > 1 show a dra- 
matic increase of the spot size due to the low energy of 
the bunches. This is illustrated in Fig. |H1 where the sim- 
ulated distribution at the detector is shown for n = 1, 



n = 2 and n = 3 at a bunch energy of 32 eV. The rms 
radius for n = 1,2,3 in the x-direction is respectively 
<j\ = 65 /im, <72 = 471 firm and 03 = 635 /im. The spot 
size due to n = 2 and n = 3 is much larger than the spot 
size due to n = 1 and even if the probabilities for n = 2, 3 
are low compared to that for n = 1 they will play a small 
role in the final spot size. The asymmetry present in the 
spots of Fig. [5] (b) and (c) is due to the fact that the ini- 
tial ionization volume is not spherical {a x < a y ). In case 
of three particles per bunch, this effect is still noticeable 
but more smeared out, see Fig. [§] (c). 
The sum distribution D(x) is given by: 

3 

D(x) = ^ n Pn(n, n)D n (x), (12) 

n=l 

where D n (x) is the distribution of the final spot for 
n = 1, 2, 3 obtained from GPT simulations and p n is the 
Poisson probability used as a weight. The sum distribu- 
tion is fitted to a 2-dimensional Gaussian. The stan- 
dard deviation a S i m of the fitted Gaussian along the 
x-direction at different V neg is also showed in Fig. [5] 
(dashed line). We find that space-charge forces have a 
very small effect on the waist-scan. The dashed line in 
Fig. [SJ drawn for an effective source temperature of 3 
mK, marginally improves the match with the experimen- 
tal data (x 2 = 1.1). While noticeable, it is uncertain if 
the resulting broadening can be actually extracted from 
the data since it is on the order of the accuracy of the 
measurement. Further investigations are required in or- 
der to quantify the influence of space charge. 

VII. CONCLUSION 

We investigated the source temperature of the UCIS, 
a new kind of ion source that can be used for FIB 
applications. Time-dependent electric fields are used to 
focus Rb + ion bunches and perform waist scans in order 
to determine the effective source temperature, found 
to be (3 ± 2) mK. The expected source temperature 
of To — 390/iK is consistent with this result given the 
error bounds; the lower value may point to residual 
distortions present in the beam line. The result also 
weakly depends on space-charge forces. 
From Eq. ((4]), the rms reduced emittance is calculated 
to be e r — 1.4 x 10~ 8 m rad y/ eV for an effective source 
temperature of 3 mK and an initial source size a Xi — 38 
/im. The only other omittance measurements of an 
ultracold ion source was performed by Hanssen et al*£. 
They measured the temperature of an ultracold Cr + 
source with a different method: observing how the size 
of an unfocused ion beam increases due to the source 
temperature when lowering the energy of the beam. 
The effective source temperature was varied by tuning 
the ionization laser to a lower wavelength, otherwise 
the source temperature would be too low to give an 
appreciable effect on the final spot size. Their measured 
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(a)n = 1 (b)n = 2 (c)n = 3 



FIG. 9. Particle distribution on the detector in GPT simulations depending on the number of ions in the bunch n. The plots 
are 2D histograms. The scale is the same in the three images: 3.2 by 3.2 mm 2 . The horizontal axis is the x direction and the 
vertical axis is the y direction. The rms radius for n = 1, 2, 3 in the x-direction is respectively <ti = 65 /im, 02 = 471 fim and 
03 = 635 ^m. The absolute number of counts per pixel is not of interest here and it is not given but the color code ranges from 
white (no count) to black (maximum number of counts). The asymmetry present in the spot size is due to the fact that the 
initial ionization volume is not symmetric. 



reduced emittance is a factor 23 smaller because they 
found a lower source temperature (Chromium has a 
slightly smaller Doppler temperature^) and used an 
initial source size a Xi — 5 /im. It is possible to create 
a smaller emittance by reducing the initial size of the 
source, but the difference in temperature can not be 
compensated. Due to the uncertainty in our measured 
source temperature, in fact the two independent mea- 
surements overlap within two standard deviations. This 
confirms that the effective temperature of the UCIS is 
indeed close to that of the laser-cooled atoms, which 
is an essential ingredient to achieve high brightness 
with the UCIS. Moreover, we presented a new method 
focusing Rb + bunches with time-dependent fields. 
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